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We compute the branching ratios for the rare top quark decays t ^ cy and t — >■ cZ mediated 
by effective flavor changing neutral couplings at the next-to-leading order in QCD including the 
effects due to operator mixing. After re-suming contributions of order [as log(A/mt)]", where A is 
the scale at which the effective operators are generated, at leading log level using RGE methods, 
we compute finite matrix element corrections and study the effects of experimental kinematic cuts 
on the extracted branching ratios. We find that the t ^ cy decay can be used to probe also the 
effective operators mediating t —> eg processes, since these can naturaly contribute 10% or more to 
the radiative decay. Conversely, any experimental signal of t ^ eg would indicate a natural lower 
bound on t ^ cZ, 7. 
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The standard model (SM) predicts highly suppressed 
flavour changing neutral current (FCNC) processes of the 
top quark (t — >■ cV, V = Z,j, g) while new physics be- 
yond the SM (NP) in many cases lifts this suppression 
(for a recent review c.f. [1]). 

Top quark FCNCs can be probed both in production 
and in decays. Presently the most stringent bound on 
the Br{t cZ) comes from a search performed by the 
CDF collaboration Br{t cZ) < 3.7% at 95% C.L. 
[5]. For the photonic decay, the most stringent bound 
was put forward by the ZEUS collaboration Br{t — > 
C7) < 0.59% at 95% C.L. [3]. On the other hand the 
most stringent present limit on Br{t — eg) comes from 
the total cross-section measurement of CDF and yields 
Br{t eg) < 0.57% at 95% C.L. g]. The LHC will be 
producing about 80, 000 events per day at the luminos- 
ity L = lO'^'^cm^^s^^ and will be able to access rare top 
decay branching ratios at the 10~^ level with 10fb~^ [5]. 

Recently the t ^ eV decays mediated by effec- 
tive FCNC couplings have been investigated at next- 
to-leading order (NLO) in QCD and it was found that 
t ^ eg receives almost 20% enhancement while correc- 
tions to the t — >■ C7, Z branching ratios are much smaller. 
However, the authors of 6J only considered a subset of 
all possible FCNC operators mediating t eV decays 
at leading order and furthermore neglected the mixing 
of the operators induced by QCD corrections. In the 
case of i — 7- C7 decay in particular the QCD corrections 
generate a nontrivial photon spectrum and the correct 
process under study is actually t — >■ ejg. Experimental 
signal selection for this mode is usually based on kine- 
matical cuts, significantly affecting the spectrum. The 
validity of theoretical estimates based on the completely 
inclusive total rate should thus be reexamined. Finally, 
renornialization effects induced by the running of the op- 
erators from the NP scale A to the top quark scale are 
potentially much larger than the finite matrix element 



corrections. Although these effects are not needed when 
bounding individual effective FCNC couplings from indi- 
vidual null measurements, they become instrumental for 
interpreting a possible positive signal and relating the 
effective description to concrete NP models. 

In this Letter we present the results for the NLO QCD 
corrections to the complete set of FCNC operators medi- 
ating t eV decays already at the leading order includ- 
ing operator mixing and renormalization effects. Finally 
in the case of i — > C7 we study the effects of experimental 
kinematical cuts on the extracted branching ratio limits 
at NLO in QCD. 

Basis of operators. In writing the effective top 
FCNC Lagrangian we follow roughly the notation of 
ref. ^[Zj. Hermitian conjugate and chirality flipped op- 
erators are implicitly contained in the Lagrangian and 
contributing to the relevant decay modes 
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= ^ E blj,OYi, + ^afOl + {L^R) + h.c., 

V=g,i,Z 

(1) 

where O^^^;,^ - gvV^,qL.,RT-a^^nR^L. OIj,^ = 
9zZf,qL,R7''tL,R, q = c{u), qR.L = (1 ± 75) g/2, c r^,. = 
«[7m,7!^]/2, gz = 2e/ sm29w, 9-, ^ e, gg = \/al47r. Fur- 
thermore V{A,Z)^, = a^K - d^V^, Gl, = d^Gl - 
d^Gl + gfabcG^Gl, and are the Gell-Mann matri- 
ces in the case of the gluon and 1 for the 7, ^. Finally 
V = 246 GeV is the electroweak condensate and A is the 
effective scale of NP. 

Note that in principle, additional, four-fermion opera- 
tors might be induced at the high scale which will also 
give (suppressed) contributions to t — ^ cV processes. On 
the other hand, such contributions can be directly con- 
strained via e.g. single top production measurements and 
we neglect their effects in the present study. 

Operator renormalization. The QCD virtual cor- 
rections to effective operators in eq. ([l]) involve ultra 
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violet (UV) divergencies. These are cancelled exactly 
in the matching procedure to the underlying NP theory. 
The remaining logarithmic dependence on the matching 
scale can be resumed using renormalization group (RG) 
methods. The RG running is governed by the anomalous 
dimensions of the operators. The operators Of ^ do not 
mix with the others under QCD renormalization. They 
have identical anomalous dimensions of 7^ ^ = agCp/T^, 
where Cp — 4/3. We assemble the remaining six oper- 
ators into two vectors Oi = {O] ,Of ,Of)'^ , i = RL.LR 
which again do not mix with each other under QCD run- 
ning. The corresponding one-loop anomalous dimension 
matrix is the same for both chiralities and reads 
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Figure 1: Effective tc^ and tcZ couplings at the top mass 
^2) scale nt ~ 200GeV induced by the RGE running of the teg 
coupling from the NP scale A. 



where Ca = 3 and = sin^ 6w — 0.231 denotes the 
square of the sinus of the Weinberg angle. Depending on 
the nature of new physics which generates these dipole 
operators at the scale A, the relevant (LR) operators 
might explicitly include a factor of the top mass (i.e. 

by re-defining operators as Olr = {mt/v)OLR) and its 
running can then be taken into account by adding 6Cf 
to the diagonal entries of 7^. As we shall demonstrate, 
this effect is numerically not important for the interest- 
ing range of couplings and scales, which can be probed 
at the Tevatron and the LHC. We are interested in par- 
ticular in the mixing of the gluonic dipole contribution 
into the photonic and Z dipole operators. For the case, 
without the top mass effects, the LR and RL operators 
receive identical corrections and the effective couplings 
at the top mass scale read 

h■l{^^t)^v'''hJ{K) + ^ (77'=^ - ^^^) fef (A) , (3a) 
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(r,«i-^«^)fef(A),(3b) 



where is the top mass scale, 77 — as(A)/as(/it), 
Ki = 4/3/3o, ^2 = 2/3/3o and /3o is the one-loop QCD 
beta function. Assuming that no new colored degrees 
of freedom appear below the UV matching scale which 
would modify the QCD beta function, it evaluates to 
/3o = 7 above the top mass scale. If we include the top 
mass running in the RGE of LR operators, then K1.2 are 
modified to ki = 16/3/3o, ^2 = 14/3/3o. 

We illustrate the effect of the RGE running in Figure 
[1] where we plot the values oi b]'^ at the top mass scale 
fit — 200 GeV induced solely by the mixing of the gluonic 
dipole contribution, produced at the UV scale A. We see 
that for NP matching scales above 2 TeV the induced 
contributions to fo^H b.l around 10% of the 
in the UV. On the other hand, due to cancelations in 
the RGE equations for the fef^^^, these receive much 
smaller corrections (below 1 % for the interesting range) . 



Including the top mass renormalization reduces the in- 
duced corrections to the coupling. However for UV 
scales of a couple of TeV or below, this effect is negligible. 

Matrix element corrections. To consistently de- 
scribe rare top decays at NLO in as one has to take 
into account finite QCD loop corrections to the matrix 
elements (97 1 Oi \t) and {qZ\ Oi \t) evaluated at the top 
mass scale as well as single gluon bremsstrahlung correc- 
tions, which cancel the associated infrared and collinear 
divergencies in the decay rates. Contributions due to the 
C'lr RL have already been computed in ref. 6 . Here we 
present results for the Of ^ as well as for the admixture 
of the gluonic dipole operators O^. The t — )■ cV rate for 
V = Z,^ with implicit operator chirality assignments (a 
and b stand for and b^ji or an and buL) reads 
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[2Re{a''*b^}Tl + 2Re{b'' * a"" }r^, 



(4) 
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where a"^ = 0. At the tree level only P^ 



b,ab 



■ I A4 

contribute 

and without imposing any kinematical cuts we derive 
them in d = 4 + e dimensions and for massless charm 
quark as 

rr'-^^ffy(l - rvrni + '-){! rv)^2(2 + e + ry) , 



167r 



16^ 



.91(1 



2rz^ 



- (2 + e)rz] , 



-m-rzYiS + e)., 



(5a) 

(5b) 
(5c) 



where rz — {raz/ratY' and = 0. We first study as 
virtual and real bremsstrahlung corrections to these con- 
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Table I: Numerical values for coefficient functions corre- 
sponding to the inputs rrit — 172.3 GeV, mz = 91.2 GeV, 
sin^ 9w 



Xb — 


2.36 


— 


1.44 




= 1.55 






Vb = 


-17.90 




-10.6 


8 Vab 


= -10.52 


Va = 


= 0.0103 


Vbg = 


-- 3.41 


Vag ~ 


2.80 


Vbg 


= 2.29 


Vag 


= 1.50 



tributions leaving operator contributions aside. We 
obtain 



An 



2 log 



(^^)+8-8Li2(ry) 
2(8 + ry) 



4 log(l - ry) \og{rv) — log(l - ry 

3 2 + ry 

Ary{2^2ry-rl) 16 - llr^/ - 17r^ 

(l-ry)2(2 + ry) "^^""^^ 3(l-ry)(2 + rv) 



(6) 




- 41og(l - rz) \og{rz) 



■ - - 1 n A 4rz(l + r2)(l-2rz) 

-2—^ log 1 - '-2 - -r, \2^^ ,0 \ log 

l + 2rz (1 - r2)^(l + 2rz) 

5 + 9rz-6r| „^ . , , 47r2 
(l ~ rz)[l + 2rz) 3 
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The as corrections to Tl in eq. Q have been derived 
before [B] while the remaining two expressions are a new 
result. We confirm the result for r^", which can most 
easily be obtained from the corresponding B — >■ Xgj ex- 
pressions, which have been known for some time [8]. On 
the other hand we find disagreement with the nonzero 
rz dependence reported in (H). We have crosschecked 
our result with the corresponding calculation done for 
a virtual photon contributing to B — Xs£^i~ In 
the case of Tf our result agrees with the correspond- 
ing calculation of i — ?> bW at NLO [10]. In fact, when 
considering the branching ratio Br{t — >■ cZ) approxi- 
mated as r{t — >■ cZ)/T{t — T' bW), as corrections cancel 
almost exactly for ^ contributions leaving corrections 
of (as/47r)6Ci.(rvy - rz)[l + C(rvi/,rz)] = -0.0006 and 
(as/47r)(CF/3)[l -hO(riv,rz)] = 0.001 respectively |14j. 
In particular, the cancelation in b^ contributions is even 
more severe than reported in ^ . 

Including contributions from gluonic dipole operators 



the analytic expressions are rather lengthy and will be 
presented elsewhere [11] . For the case of the Z boson we 
find for the completely inclusive rate 
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Re{b^*a^}(xai>+^CFyai 



A4 



+2-Re{b^*b3}^Cpy,g - 2lTn{b^*ba}^CFy,g 
47r 47r 



where the numerical Xi,yi, jji coefficients are given in Ta- 
ble |l] The resulting Ug corrections to Br{t cZ) are 
at the order of a few percent for ~ \b^\,\a^\. In 
the case of t — 07 at NLO the process in question is 
t — >■ cyg and involves three (one almost) massless par- 
ticles in the final state. Virtual matrix element correc- 
tions contribute only at the soft gluon endpoint {Eg = 0) 
and result in non- vanishing b'^b^ interference contribu- 
tions. They involve IR divergencies which are cancelled 
by the real gluon emission contributions. These also pro- 
duce non- vanishing |6^p contributions, and create a non- 
trivial photon spectrum involving both soft and collinear 
divergences. The later appear whenever a photon or a 
gluon is emitted collinear to the light charm jet. In the 
analogous B Xgj decay measured at the B-factories 
the photon energy in the B meson frame can be recon- 
structed and a hard cut (E'™*) on it removes the soft pho- 
ton divergence. The cut also ensures that the B Xgg 
process contributing at the end-point E'.^ = is sup- 
pressed. On the other hand, in present calculations the 
collinear divergencies are simply regulated by a non-zero 
strange quark mass, resulting in moderate log(r?T,s/mf,) 
contributions to the rate. The situation at the Tevatron 
and the LHC is considerably different. The initial top 
quark boost is not known and the reconstruction of the 
decay is based on triggering on isolated hard photons 
with a very loose cut on the photon energy (a typical 
value being > 10 GeV in the lab frame [Hj). Isola- 
tion criteria are usually specified in terms of a jet veto 
cone AR = ^/Arp~+~K(j)^ where Arj is the difference in 
pseudorapidity and A(j) the difference in azymuthal angle 
between the photon and nearest charged track. Typical 
values are AR > (0.2 - 0.4) [13]. We model the non- 
trivial cut in the top quark frame with a cut on the pro- 
jection of the photon direction onto the direction of any 
of the two jets {6r = 1 — p-,, • pj/E^Ej). The effects of the 
different cuts on the Dalitz plot are shown in Figure [2] 
Since at this order there are no photon collinear divergen- 
cies associated with the gluon jet, the Sr cut around the 
gluon jet has numerically negligible effect on the rates. 
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Figure 2: The t — >■ c-^g Dalitz plot. Contours of constant 
photon and gluon infrared and collinear divergent contribu- 
tions are drawn in red and blue respectively. The collinear 
divergencies appear at the horizontal and vertical boundaries 
of the phase-space, while the IR divergencies sit in the top 
and right corners. The cuts on the photon energy correspond 
to vertical lines, the cuts on the gluonic jet energy to hori- 
zontal lines. Green lines correspond to cuts on jet veto cone 
around the photon. 
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Figure 3: Relative size of Os corrections to the Br{t — >■ C7) 
at representative ranges of 5r and Ej^f Contours of constant 
correction values are plotted for fe^ = (gray, dotted), 6® = h'' 
(red) and V — —b"' (blue, dashed). 



On the other hand the corresponding cut on the charm 
jet - photon separation does not completely remove the 
collinear divergencies in the spectrum. However, they 
become integrable. The combined effect is that the con- 
tribution due to the gluonic dipole operator can be much 
more pronounced than in the case of i? — ?> X^^. 

The analytic formulae for the T^^ ^ and F^^ with full 
6r and i?™* dependence are rather lengthy and will be 
presented elsewhere [TT] . In Figure [s] we show the 6^ 
induced correction to the tree-level Br{t 07) for rep- 
resentative ranges of 5r and Ej^^. We observe, that the 



contribution of W can be of the order of 10 — 15% of 
the total measured rate, depending on the relative sizes 
and phases of O^'^ and on the particular experimen- 
tal cuts employed. Consequently, a bound on Br{t — >■ 07) 
can, depending on the experimental cuts, probe both h^^^ 
couplings. 

In summary, QCD corrections to FCNC coupling me- 
diated rare top decays can induce sizable mixing of the 
relevant operators, both through their renormalization 
scale running as well in the form of finite matrix ele- 
ment corrections. These effects are found to be relatively 
small for t ^ cZ decays. On the other hand the accurate 
interpretation of experimental bounds on radiative top 
processes in terms of effective FCNC operators requires 
the knowledge of the experimental cuts involved and can 
be used to probe O^j^ contributions indirectly. 
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